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Abstract 

Integral identities that hold between "desired" and "comparison" solutions of the radial Dirac 
equations for scattering precesses are considered. Applications of these identities are discussed, 
particularly the determination of bounds to variational calculations of iT-matrix elements. 

Relativistic effects in atomic scattering, and indeed in many quantum scattering processes, 
can usually be treated as small corrections, and handled by perturbation theory with respect to 
the non-relativistic (Schrodinger) results. However, for some processes, such as the scattering of 
electrons by atoms and molecules, and for many nuclear scattering processes, it is often convenient 
or even necessary to use the Dirac equation directly. This is because the relativistic kinematics and 
spin effects are then automatically "built in" (see, for example, ref. [1 - 3] and citations therein). 
In addition to this, relativistic corrections to the dynamics (i.e. corrections to the static potential) 
may need to be taken into account. 

In the Dirac formalism, the description of the scattering of a fermion (such as an electron or 
positron) by a target (such as a neutral atom) is often reducible to the solution of the radial Dirac 
equations (h = c = 1) 

f{r) + !f(r) = ( E + m- V(r))g(r) + X(r), (1) 
r 

9'{r) - -g(r) = -(E-m- V(r))f(r) + Y(r), (2) 
r 

where f(r) and g(r) are the usual reduced radial coefficients of the "large" and "small" components 
of the Dirac spinor for the incident fermion of mass m and energy E. The terms X(r) and Y(r) 
include exchange effects as may be applicable. We shall restrict our discussion to the class of 
potentials V(r), and corresponding exchange terms X(r),Y(r), such that the solutions f(r),g(r) 
satisfy the following boundary conditions: 

Ur = 0) = g K (r = 0) = 0, (3) 



f K (r -» oo) ~ ^(fc) sin (kr - + B K (k) cos (At - , (4) 

Ur oo) ~ ir^ = O 4 ^ cos ( fcr " 4) " ^ sin ( fer ~ 4)) ' (5) 

where k 2 = E 2 - m 2 , and k = j + \ if j = I — \ ("spin down") and k = — (j + ±) if j = £ + ± 
("spin up"). This means that the potentials are short range and not overly singular at the origin. 
Specifically, a sufficient condition is that the limit of r 2 V(r) be zero as r — > oo and as r — > 0. If the 
potential V(r) is long-range, that is contains a Coulombic contribution, then the sine and cosine 
functions in (4) and (5) would be replaced by the corresponding Coulomb functions. 
The asymptotic forms (4) and (5) can be written in the equivalent form 

f K {r -> oo) ~ C K (k) sin (kr -t^ + r/ K (fc)) , (6) 
g K (r -> oo) ~ — C K (fe)cos (At + Tj K {k) j , (7) 

+ 771 V Z / 

where ry K (/c) are the scattering phase shifts, while 

A K (k) = C K (k) cos rj K ( k) and B K (k) = C K (k) smrj K (k). (8) 

The asymptotic normalization constants C K (k) (or, equivalently, the constants A K (k),B K (k)) 
may be chosen to be anything that is convenient. Some common choices are C K (k) = l,C K (k) = 
sec r] K (k), etc. The scattering cross sections or polarization parameters are then calculated from 
the phase shifts r) K (k) [1,2]. 

In non-relativistic (Schrodinger) scattering theory perturbative effects can be taken into ac- 
count by using the integral identity between a "given" and "comparison" solution first obtained by 
Hulthen [4] and later elaborated by Kato [5] and others. This integral identity can also serve as the 
basis for approximate variational solutions to the scattering equations [4-7], and for determining 
bounds on approximate calculations of scattering parameters [8,9]. 

Evidently, analogous results can be written down in the Dirac formalism of scattering theory, 
as we now proceed to discuss. Thus, suppose f(r), g(r) are solutions of a "trial" or "comparison" 
problem, corresponding to X(r), Y(r) and the potential V(r), namely 

l\r) + -7(r) = {E + m — V{r))g{r) + X(r), (9) 
r 

g'(r) - ^g{r) = -(E-m- V(r))J(r) + Y(r), (10) 

with 

J K (r - oo) ~ C K (k) sin (kr - l\ + rj K (k)) , (11) 
2 



g K (r^oo)~^—C K (k)coB(kr-e^ + Ti K (k)^ , (12) 
Straightforward manipulations of the equations (1), (2) and (9), (10) result in the identity 

^ {1 9 - ffl = (V-V) (// + gg) + Xg-Xg + jY- f¥. (13) 

Integration of Eq. (13) leads to the result 

[f(r)g(r) - f(r)g(r)] * = f dr [(V - V)(ff + gg) + Xg-Xg + Yj- Yf] , (14) 

J o 

where f = f(r), etc. in the integrand of Eq. (14). If we now make the replacements / = /+(/ — /), 
etc., Eq. (14) can be rewritten in the form 

A = [ R dr Uv - V){f + g 2 )} + [* dr(V - V) [(/ - 7)7 + G? - g)g] 
Jo 1 J Jo 

+ I* dr [(X - X)g - (Y - Y)J] + f dr \X(g -g)- Y(f - J)] , (15) 
Jo Jo 



where 



A = [f(r)g(r) - f(r)g(r)]* ~ -J— CCsin(ry K - r, K ) = -^—{BA - AB), (16) 
L u h + m h + m 



and the symbol ~ indicates that R has been taken to be sufficiently large that the asymptotic forms 

(6), (7), (11) and (12) apply (we can take R — > oo). 

The integral identities (14) and (15) relate the phase shifts r] K (or, more generally, functions 

of these, such as the K-matrix elements, K K = tan?? K , T-matrix elements, T K = e irjK sin?7 K , etc.) 

to the "comparison" phase shifts fj K (or corresponding functions thereof). This is clear from the 

explicit form of A for given choice of asymptotic normalization, that is, choice of A and B or 

_ 

alternatively C. For example if C = C = 1 then A = — sin(rj„ — r] K ), or if A = A = 1, B = 

E + m 



tanr/ K = K K , B = tanr/ K = K K then A = — (K K - K K ), etc. 

hi + m 

The integral identities (14) or (15) can be used for various purposes, some of which we discuss 
briefly in what follows: 

1. Formal results. 

If we take V = 0, and the corresponding free incident wave solutions of Eqs. (9) and (10) for 
f K and g K , then the identity (14) (with A = A = 1, B = tanr/ K = K K , B = 0) gives the well-known 
integral expression for the K— matrix elements, 



K K = - 



E + m ' /> ' 



k 



[ dr[V(ff + gg)-Xg + Yf]. (17) 
Jo 



This is often used for extracting the phase shifts from numerical solutions of Eqs. (1) and (2). 



2. Perturbative calculations. 

A not-infrequent situation is that the potential V can be written in the form V = V + Vi, 
where Vq is a dominant (and/or easily solvable) interaction term (such as the electrostatic potential 
in atomic scattering), and V\ is a small "correction" term. Then, obviously, if V = Vq, f and g are 
known (or easily obtainable), while V — V = V\ can be handled perturbatively. Taking f = f and 
g = g in lowest order on the right-hand-side of Eq. (15), one can use that equation to evaluate r\ K 
in terms of rj K plus a lowest order perturbative correction (which is given by the RHS of Eq. (15) 
with f = f and g = ~g). In general the perturbation may be in V only, or in X and Y, or both (see, 
for example, ref. [10]). 

3. Variational approximations. 

In some instances it may be useful or necessary to approximate the solutions of (1) and (2) 
variationally. For example, one may wish to have analytic representations of the solutions (recall 
that, with rare exceptions, Eqs. (1) and (2) are not analytically solvable). In such cases, one can use 
a variational approach, in which the desired (unknown) solutions /(r), g{r) are approximated by 
analytic trial forms /(r), g(r) that contain adjustable parameters ctj (j = l,...,n p ). The identity 
(15) can be used to choose these parameters ctj in a variationally optimal way. We illustrate this on 
the case X = Y = X = Y = 0, and normalization choice C K = C K = 1, in which case the identity 
(15) can be written as 

■sm(r] K -r ]K )=I[f,g]+n 2 [f,gJ,g], (18) 



E + m 

where 

2 ) (!»") 

r-R 



I[f,9} = [ R dr(V-V)(f 2 + g 2 
Jo 



[ dr 
Jo 



fU±-^g + (E- m -V)A-gUf r + ^f-(E + m -V)g 



,(1%) 



(20) 



and where we have used the identities 

{jr + 7) 7(0 " {E + m " V{r)) W) = (F(r) " F(r)) W) ' 

" 7) v{r) + ( E ~ m - y ( r )) 7( r ) = - ( y ( r ) - F ( r )) (21) 

in rewriting (19a) in the form (19b). The term 1Z 2 is a "remainder" that is given by the expression 
K2 [/, 9, 7, g] = j R dr (V - V) [(/ - 7)7 + (9 ~ M , (22) 

JO 
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which is second order in the "small" quantities / — f,g — g, and V — V. Usually we take R — > oo 
in these integrals, and this shall be done in the rest of this paper. 
From Eq. (18), if we neglect IZ2, it follows that 

r, K (k) ~ fj K (k) - sin" 1 1[f,g] = r,^(k, aj ), (23) 

where r]i App '\k,aj) is the approximate value of rj K (k) for any given k and k. Note that an explicit 
knowledge of the comparison potential V is not necessary to evaluate rj^ pp '\k, ctj), that is, it is 
only necessary to choose the trial functions /, ~g. This is evident from Eq. (19b), in which, as can 
be seen, V does not appear explicitly. Of course, we want to choose the adjustable parameters ctj 
of f(r, ctj) and ~g(r, ctj) in such a way that ^ App0 (ib,a 3 is as close to rj K (k) as possible. In other 



words, we wish to minimize 
d 



n K {k)-n {App -\k,a 3 



with respect to aj. Since 



rj K (k) - ■qi App ) (k,a j 







r, K (k)- V ^(k,ctj) =- K — L _„(AppO (A . >a)j (24 ) 



V K (k) - rf App -\k,a ji 



we see that a condition for a minimum of \rj K (k) — r/i App '\k,aj)\ is that 

d 



daj 



4 A ^-\k,a J ) = 0, (25) 



as happens also in the corresponding Schrodinger theory. The resulting optimal values, a° pt , of the 

adjustable parameters aj (which include the trial value rj K of the phase shift, or of K K = tan 77,, if 

the normalization A = 1,B = tanr/ K is used, etc.), are then substituted into Eq. (23) to yield the 

optimal variational approximation, r/i App '- 1 (k, a° pt )> to ij K (k) (or K K = tan(r/ K ), etc.), corresponding 

to a minimum value of 1 72-2 1- (Strictly speaking, minimuma of \rj K (k) — rf App '\k,aj)\ may occur 

d 

at points in parameter space where — — r]^, App '\k, aj) is undefined (i.e. cusps rather than smooth 
minima) or at boundary points of the domain of parameter space. Such possibilities must be kept 
in mind and investigated, if necessary.) 

4. Bounds on scattering parameters. 

In approximate calculations of scattering parameters (phase shifts, K-matrix elements, etc.) 
neither the sign nor the magnitude of the difference between the (unknown) exact and approximate 
value is known. However, for the case X = X = Y = Y = 0, if we write V = V + SV, f = f + 5f 
and g = g + 5g, where 5V, 8f, 5g — > 0, then (with the choice of asymptotic normalization C = 1) 
Eq. (14) implies that 

pr-(^-%) = - dr5V(f + f), (26) 
h + m J 



where we have kept only the first order terms in infinitesimal quantities, and so set s'm(5r]) = 5rj. 
Equation (26) shows that if V — > V from below, i.e. if SV = V — V < 0, then fj K > r] K and 
vice- versa (as happens also in Schrodinger theory). This property can be used to set up a scheme 
in which approximate calculations of phase shifts approach the (unknown) exact values from above 
(or below), provided that the trial solutions are so chosen that the corresponding trial potentials 
approach the exact one from below (or above). 

Although, as already stated, in general it is not possible to evaluate either the sign or the 
magnitude of the remainder term 1Z 2 , Eq. (22), it is possible, in some cases, to determine calculable 
bounds B on 1Z 2 of the form 

\n 2 [f,gJ,g]\<B[VJ,g}. (27) 

This, together with Eq. (18) (or its equivalent with other asymptotic normalizations), leads to upper 
and lower bounds on the scattering parameters. We illustrate this on the potential scattering case 
(X = Y = X = Y = 0), and the choice of asymptotic normalization A = A = 1, B = tanr/ K = 
K K , B = tanr/ K = K K , whereupon Eq. (14) becomes 

(K K -K K ) = I[f,g}+TZ 2 , (28) 



E + m 

where I[f,g] is given in Eq. (19). 

We write the remainder term, Eq. (22) in the form 

n 2 = n 2L + n 2S , (29) 

where 

U 2L = / drAVF{r)J(r), K 2S = / drAVG(r)g{r), (30) 
Jo Jo 

with A V = V—V, F = f—f&ndG = g—g. Then, using the Schwartz inequality (s,t) 2 < (s,s)(t,t), 
it follows from Eq. (29) that 



where 



\R-2l\ < ciFbj, \R- 2 s\<<iGbg and so \1Z 2 \ < a F bj + a G bg, (31) 

poo />oo 

b 2 J = I drp-\r)\AV(r)J(r)\\ 6| = J dr p~ l (r)|AV (r)g(r)\ 2 , (32) 

a%= drp(r)\F(r)\ 2 , a 2 G = drp(r)\G(r)\ 2 , (33) 
Jo Jo 

and p(r) is an arbitrary, positive weight function (but such that all the indicated integrals exist). 
For example, p(r) might be chosen to be |AV(r)|, or some other positive function, possibly with 

6 



adjustable parameters, such that the indicated integrals exist. For a given choice of trial functions 
/, g, the expressions bj and bg of Eq. (32) are calculable (remember that V need not be known 
explicitly, in light of the identities (20) and (21)). 

It remains now to determine bounds on ap and aa (which are not calculable since F = f—f and 
G = g — g are not known). One way that such bounds can be obtained is from the integral equations 
for the radial Dirac functions (written here for the present choice of asymptotic normalization 
A = A = 1 so that B = K = tan 77 and B = K = tan rj) : 

l-CO 

/«(r) = ui(r)+ / dT>U{T>)[G\\rS)Ur>) + G?{rS)g K {T>)], (34) 
J o 



g K (r)=u 2 (r)+ / dr'U{r') [Gf (r, r')f K (r') + Gf (r, r')g K {r')] , (35) 
Jo 

where the Green functions Gf 3 are defined by 

Gf (r, r') = -v a (kr)u b (kr') r 1 < r 

= \u a {kr)v b (kr') r' > r (36) 

and U(r) = (E + m)V(r). The functions u a , v a are defined in terms of the usual Ricatti-Bessel and 
Ricatti-Neumann functions [11], 

je(kr) = krje(kr) ~ sin(/cr — £—) and h(,(kr) = krni{kr) ~ — cos(/cr — ^)-, (37) 

namely 

ui(kr) = jg(kr) vi(kr) = hi(kr), (38) 

k k 

u 2 (kr) = cr K — j-Akr) v 2 (kr) = o K — n-(kr), (39) 

hi + m h + m 

where £ = £ — a K and a K = - — - is the sign of k. Similar integral equations can be written down for 

\k\ 

f and g, hence also for F = f — f and G = g — g, specifically 

POO 

F K {r) = ~F K (r) + / dr'U(r') [G\\r, r')F K {r') + Gl 2 (r,r')G K (r')] , (40) 
J o 

/>oo 

G K (r) = G K {r)+ / dr'U(r') [Gf (r,r')F K (r') + Gf (r,r')G K (r')] , (41) 
Jo 

F K (r)= / dr'(U(r>)-U(r>))[Gf(r,r>)J K (r>) + Gf(r,r')g K (r>)}, (42) 

f'OO 

G K (r)= / dr'{U(r')-U(r'))[Gf( r y)J K (r') + Gf(ry)g K (r')]. (43) 
Jo 



where 



Note that F K (r), G K (r) are known functions, for given trial functions f K and g K , since G" 6 (r, r') 
are known. We stress that the explicit form of the trial potential U(r) need not be known in Eqs. 
(42) and (43) because of the identities (20) and (21). Thus, only the trial functions f K (r,ctj) and 
~g K (r,ctj) need be specified. 

Now, multiplying Eq. (42) by p(r)F*(r), integrating over r, and making repeated use of the 
Schwartz inequality gives the result 

a F <a T + a F g n + a G g 12 , (44) 

and similarly 

a G < olq + a F g 2 i + a G g 22 , (45) 

where a-p and a^, defined as in Eq. (33), are calculable since F and G are known. The factors gij 
are given by 

g f. = / / drdr , p(r)\G i j(r,r')U(r')\ 2 p- 1 (r'). (46) 
Jo Jo 

The generalized Schwartz inequality 

J J drdr's(r)Q(r,r')t(r') < j dr\s{t)\ 2 j dr'\t(r')\ 2 j J drdr'\Q(r,r')\ 2 , (47) 

was used in obtaining the results (44)- (46). 
From Eqs. (44) and (45), it follows that 

a F < p [(! - 522)0^ + 5i20g] = ( 48 ) 



<i g < ^ [(1 - 911)0^ + 521 a ¥ ] = B G , (49) 



provided that 

gu<l (* = 1,2) and 23 = (1 - 5 n)(l - 522) - £u 2 £T2i > 0. (50) 

From the definition (46) of g^, and that of the Green functions (36), it is clear that the conditions 
(50) are, for given k (i.e. given energy of incidence), restrictions on the strength of the potential 
V(r). That is, the potential must be sufficiently weak for the inequalities (50) to be met. Note, 
however, that since the Green functions contain the factor — , g^ will generally decrease with 
increasing k. This means that a given potential V(r) may be such that the inequalities (50) might 
not hold when k is small (low-energy scattering) but will hold for higher values of k. 

Replacing the expressions a F and a G by their bounds (48) and (49) in Eq. (31) then leads to 
the inequality 

\TZ 2 \<B F bj + B G bg = B 2 , (51) 



and hence to the following simultaneous upper and lower bounds on the (unknown) exact if-matrix 
element: 

x ( App0 _ E + m^ < ^( Exact) = tan ^ < K ( App . } + E+H^ (52) 

where 

K^=K-^±^I[f,g\. (53) 

The definition (53) is the analogue of that of Eq. (23) for the present choice of asymptotic nor- 
malization A = A = 1. Note, again, that the bounds (52) hold provided that all the integrals that 
enter into the expressions for K^ App ^ and B2 exist, and that the inequalities (50) apply. 

We stress that the bound £>2 of Eq. (51) (with Eqs. (32), (48) and (49)) is expressible in 
terms of V(r) and the trial functions f K (r,otj) and g K (r,ctj), hence it is ultimately a function of 
the adjustable parameters, that is 82(0 j). These parameters may be chosen in accordance with 
the variational prescription (25), or such that the upper and lower bounds are as close as possible, 
i.e. such that B 2 (aj) is a minimum. These two prescriptions are not the same but, for sufficiently 
flexible trial functions, they will yield similar results. In practice, the prescription (25) is simpler to 
implement. In either case £>2 can be made as small as desirable (in the domain where the conditions 
(50) hold), provided that / and g are sufficiently flexible. 

To summarize, we have presented integral identities that hold between given and a comparison 
(or "trial") solutions for scattering calculations in the Dirac formalism. Various applications of these 
integral identities have been discussed, including their use in approximate, variational solutions of 
the scattering parameters (phase shifts or functions thereof). In particular, we have used these 
integral identities to establish rigorous and calculable bounds on the difference between the exact 
and approximate i^T-matrix elements for a wide class of potentials. These bounds can be made as 
tight as necessary, if sufficiently flexible trial functions are used. 
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